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Abstract 

We investigate the effects of hydrogen plasma treatment (HPT) on the properties of silicon quantum dot 
superlattice films. Hydrogen introduced in the films efficiently passivates silicon and carbon dangling bonds at a 
treatment temperature of approximately 400°C The total dangling bond density decreases from 1.1 x 10 19 cm -3 to 
3.7 x 10 17 cm" 3 , which is comparable to the defect density of typical hydrogenated amorphous silicon carbide films. 
A damaged layer is found to form on the surface by HPT; this layer can be easily removed by reactive ion etching. 
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Background 

Solar cells that use nanomaterials have attracted interest 
for their potential as ultra-high efficiency solar cells [1]. 
The conversion efficiency limit of a single-junction solar 
cell strongly depends on the band gap of the absorber 
layer, which is known as the Shockley-Queisser limit [2]. 
To overcome the efficiency limit, various types of 
quantum dot solar cells, such as quantum size effect 
type, intermediate band type, and multiexciton gener- 
ation type, have been proposed [3-5]. The quantum size 
effect type utilizes the phenomenon that the band gap 
of a material can be tuned by controlling the diameter 
of quantum dots, including the periodically arranged 
narrow-gap quantum dots in a wide-gap dielectric 
matrix. The fabrication of an amorphous silicon dioxide 
(a-Si0 2 ) matrix including size-controlled silicon quantum 
dots (Si-QDs) was reported by Zacharias et al. [6]. The 
size-controlled Si-QDs can be formed by annealing a 
superlattice with silicon-rich silicon oxide layers and 
stoichiometric silicon oxide layers, which is called a silicon 
quantum dot superlattice structure (Si-QDSL). Since this 
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report was published, silicon quantum dots embedded in 
various wide-gap materials, such as amorphous silicon 
carbide (a-SiC), amorphous silicon nitride (a-Si 3 N 4 ), and 
hybrid matrices, have been reported [4,7-11]. Further, the 
quantum size effect can be observed from the measure- 
ment of photoluminescence spectra or absorption coeffi- 
cients [12-14]. The Bloch carrier mobility in a Si-QDSL 
with an a-SiC matrix is higher than that in a Si-QDSL with 
an a-Si0 2 or an a-Si 3 N 4 matrix [15]. The barrier height be- 
tween a-SiC and Si quantum dots is lower than those of the 
other two materials, resulting in the easy formation of mini- 
bands [16]. Moreover, the crystallization temperature of 
a-SiC is lower than those of the other materials. Therefore, 
in this study, we focus on a Si-QDSL with an a-SiC matrix. 

High-temperature annealing above 900°C is needed to 
fabricate a Si-QDSL with an a-SiC matrix. Several prob- 
lems, such as dehydrogenation, formation of leakage 
paths by the crystallization of the a-SiC phase, and dop- 
ant diffusion in the intrinsic layer of the solar cell struc- 
ture, occur during high-temperature annealing, resulting 
in the degradation of solar cell performance. The de- 
hydrogenation problem has been addressed by hydrogen 
plasma treatment (HPT) [17]. The crystallization of the 
a-SiC phase can be prevented by incorporating a small 
amount of oxygen in the a-SiC matrix [16]. Niobium- 



Spri 



ringer 



© 2014 Yamada et al.; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons 
Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction 
in any medium, provided the original work is properly cited. 



Yamada et al. Nanoscale Research Letters 2014, 9:72 
http://www.nanoscalereslett.eom/content/9/1/72 



Page 2 of 8 



doped titanium dioxide (Ti0 2 :Nb) can be used as a 
phosphorus (dopant) diffusion barrier layer for the 
Si-QDSL solar cell [18]. Using these techniques, an effi- 
ciency of 039% has been achieved in Si-QDSL solar cells 
fabricated on insulator substrates [19]. Some researchers 
have reported the electrical properties of silicon quantum 
dot solar cells [20,21]. However, clear evidence of the con- 
tribution from Si-QDs has not yet been reported because 
of poor device quality. To improve device quality, the 
collection efficiency of the photogenerated carrier should 
be improved. For this purpose, further reduction of the 
defect density in the Si-QDSL layers and improvement 
of the p/i interface is significantly important. 

In this study, the dependence of hydrogen concentra- 
tion and defect density in Si-QDSL films on the process 
temperature of HPT was investigated. Diffusion coeffi- 
cients of hydrogen in Si-QDSLs for several treatment 
temperatures were estimated by secondary ion mass spec- 
trometry (SIMS). Hydrogen incorporation was also inves- 
tigated by Raman scattering spectroscopy. In addition, 
spin densities were measured by electron spin resonance 
(ESR) spectroscopy, and the optimal temperature was 
explored. The influence of HPT on the surface of 
Si-QDSLs was also investigated. The surface morphologies 
of Si-QDSLs after HPT were measured by atomic force 
microscopy (AFM), and the thicknesses of the surface 
damaged layers were estimated by spectroscopic ellipsome- 
try and cross-sectional transmission electron microscopy 



(TEM). The etching of the surface damaged layer was 
performed by reactive ion etching (RIE) using a tetrafluoro- 
methane and oxygen (CF 4 + 0 2 ) gas mixture. 

Methods 

Forty-period hydrogenated amorphous silicon oxycarbide 
with a silicon-rich composition (a-Sio.56Co.32Oo.12-H) /hydro- 
genated amorphous silicon oxycarbide (a-Sio.40Co.35Oo.25-H) 
superlattice was deposited on quartz substrates using 
very-high frequency plasma-enhanced chemical vapor de- 
position. The source gases were silane (SiH 4 ), monomethyl- 
silane (MMS), hydrogen (H 2 ), and carbon dioxide (C0 2 ). 
The flow rates of MMS, H 2 , and C0 2 were fixed as 1.7, 
47.5, and 0.4 seem, respectively. SiH 4 was intermittently 
flowed during the deposition of silicon-rich layers. Plasma 
power density, plasma frequency, deposition temperature, 
deposition pressure, and electrode distance were 13 mW/ 
cm 2 , 60 MHz, 193°C, 20 Pa, and 3 cm, respectively. The 
thicknesses of silicon-rich layers and stoichiometric layers 
were 5 and 2 nm, respectively. The films were thermally 
annealed at 900°C for 30 min under a forming gas (3% H 2 
+ 97% N 2 ) atmosphere to form Si-QDs. Film thicknesses of 
post-annealed samples were 250 ± 10 nm. 

After annealing, the samples were exposed to hydro- 
gen plasma to terminate dangling bond defects accom- 
panying hydrogen atoms in the Si-QDSL. The flow rate 
of H 2 , plasma power density, plasma frequency, process 
pressure, and electrode distance were 200 seem, 2.60 
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Figure 1 Depth profiles of hydrogen concentrations, (a) At 300°C for 20 min. (b) At 400°C for 10 min. (c) At 500°C for 3 min. (d) At 600°C 
for 1 min. 
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Figure 2 Arrhenius plot of diffusion coefficient of hydrogen in Si-QDSLs. The calculated diffusion coefficients of single-crystal silicon by van 
Wieringen et al. [22] and the estimated diffusion coefficients of an a-SiC thin film with hydrogen concentration of 0.4 ± 0.1 at.% by Schmidt et al. 
[23] are also described. 



W/cm 2 , 60 MHz, 600 Pa, and 3 cm, respectively. The 
treatment temperature was varied from 200°C to 600°C 
To evaluate the hydrogen diffusion coefficient in the 
Si-QDSL, the samples were treated at 300°C for 20 min, 
400°C for 10 min, 500°C for 3 min, and 600°C for 1 min. 
The depth profiles of the hydrogen concentration were 
measured by SIMS. In the measurements, Ce + ions were 
used to measure the hydrogen depth profiles. Also, the 
depth was calibrated by the etching rate of the Si-QDSL. 
Crystalline silicon was used as the standard sample to 
evaluate the hydrogen concentration. The accuracy of 
the hydrogen concentration by the SIMS measurement 
was ± 40%. In addition, for measurements of Raman 
scattering spectra and ESR, treatment temperature was 
varied from 200°C to 600°C and the treatment time was 
fixed at 60 min. 

The thicknesses of surface damaged layers formed by 
60-min HPT were estimated by spectroscopic ellipsome- 
try and cross-sectional TEM. The surface morphologies 
of Si-QDSLs after a 60-min HPT were measured by AFM. 
The etching of the surface damaged layer was performed 
by RIE using CF 4 + 0 2 gas (4% 0 2 + 96% CF 4 ). The gas 
flow rate, process pressure, and plasma power density 
were 10 seem, 4 Pa, and 0.221 W/cm 2 , respectively. The 



surface morphologies after etching were evaluated by 
AFM and spectroscopic ellipsometry. 

Results and discussion 

An average hydrogen concentration of 8.2 x 10 22 cm" 3 
was almost uniformly incorporated in the superlattice 
films before thermal annealing. After annealing at 900°C, 
the average hydrogen concentration decreased to 1.4 x 
10 20 cm" 3 . After HPT, the hydrogen concentration 
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Figure 3 Spin densities of Si-QDSLs after a 60-min HPT. 
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Figure 4 Raman spectra of Si-QDSLs after a 60-min HPT. 
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increased. Figure 1 shows the depth profiles of hydrogen 
concentrations of Si-QDSL samples treated at 300°C for 
20 min, 400°C for 10 min, 500°C for 3 min, and 600°C 
for 1 min. The oscillations with small amplitudes in the 
depth profiles are due to the matrix effect caused by car- 
bon in the Si-QDSLs. The influence of the matrix effect 
can be negligible. In addition, structure of the Si-QDSL 
is almost uniform in the depth direction. Therefore, one 
can believe the shape of the hydrogen depth profile, 
which is important to determine the hydrogen diffusion 
coefficient. The diffusion coefficients can be estimated 
from these depth profiles. The hydrogen diffusion 
process follows the diffusion equation 



dC 

dt'' 



D 



PC 

dx 2 ' 



(i) 



where D is the diffusion coefficient and C is the hydro- 
gen concentration at depth x and time t If the initial 
hydrogen concentration at all x can be assumed to zero 
except for x = 0, and if the hydrogen concentration at 



x = 0 can be fixed as C 0 for t > 0, the following solution 
is obtained: 



C(x 1 t) = Co erfc 



x 



2VDt 



(2) 



where the function erfc(z) is a complementary error 
function. The average diffusion coefficients were esti- 
mated by fitting the depth profiles with Equation 2. Red 
lines in Figure 1 indicate the fitting curves based on 
Equation 2. The calculated diffusion coefficients for each 
temperature were described by dots in Figure 2. The 
diffusion coefficient obeys Arrhenius law: 



D = D 0 exp 



AE 



(3) 



where D 0 denotes the preexponential factor, AE is the 
activation energy, and k B is the Boltzmann constant. 
From the result of the fitting by least squares method, 
D 0 and AE were estimated as 3.93 x 10" 7 cm 2 /s and 
0.81 eV, respectively. The calculated diffusion coeffi- 
cients of single- crystal silicon by van Wieringen et al. 
[22] and the estimated diffusion coefficients of an a-SiC 
thin film with hydrogen concentration of 0.4 ±0.1 at.% 
by Schmidt et al. [23] are also described in Figure 2. D 0 
and AE for single-crystal silicon and the a-SiC thin film 
are 9.67 x 10 3 cm 2 /s and 0.48 eV and 0.71 cm 2 /s and 
3.2 eV, respectively. Compared with these AE values, AE 
for Si-QDSL is relatively close to the AE for single- 
crystal Si. Such small AE indicates that the interstitial 
diffusion in Si-QDs is dominant because the thickness of 
the a-SiCO layers is too thin to work as barriers against 
hydrogen diffusion; this is due to the wide band gap and 
polar bonds of a-SiC [24]. 

From the depth profiles of Si-QDSLs for a treatment 
temperature of 600°C, hydrogen concentration was found 
to drastically decrease. Saturation hydrogen concentration 
after sufficient treatment was estimated at approximately 
1.0 x 10 21 cnT 3 , indicating that the hydrogen concentra- 
tion at the surface drastically decreases because the loss 
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Figure 5 Schematic of the structure of Si-QDSLs after HPT for the parameter fitting of spectroscopic ellipsometry. 
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of adsorbed hydrogen atoms is dominant at high 
temperatures. 

The defect densities of Si-QDSLs after 60-min HPT for 
several treatment temperatures were measured by ESR. 
The defect densities originating from silicon dangling 
bonds (Si-DBs) and carbon dangling bonds (C-DBs) were 
also estimated. The waveform separation of the obtained 
differentiated waves originating from both Si-DBs and 
C-DBs were so difficult that the ratios between the dens- 
ities of Si-DBs and C-DBs were estimated by the following 
equations [25]: 



gfa y) = Yl 



(l-x)gsi( n )+ x Ygc( n ) r x n (l-x) 12 - n 
(!-#)+ ay 12 n 1 ; 



A/si-DB — z . A/"xotal-DB 

1 + y 



and 



N c - 



DB 



i + r 



^Total-DB, 



(4) 
(5) 

(6) 



where Af Total _ DB , A^ Si _ DB , and Af c _ DB are the densities of 
total dangling bonds (Total-DBs), Si-DBs, and C-DBs, 
respectively, y is the ratio of N C -db to A^ Si _ DB and x is the 
composition ratio of C to Si. g Si {n) is the average g value 
of Si-DB surrounded by (12 - n) Si atoms and n C 
atoms, g c (n) is the average g value of a C-DB sur- 
rounded by (12-n) Si atoms and n C atoms, and g(x, y) 
is the measured g value. Although this expression is de- 
rived for an a-Sii^Q alloy system, it is believed to be valid 
for Si-QDSL with an a-SiC matrix, which can be consid- 
ered as an approximately homogeneous material, since the 
dangling bond defect density in Si-QDs is much lower 
than that of the a-SiC matrix, and the dangling bonds on 
Si-QD surfaces are passivated by the a-SiC matrix. An 
average composition ratio of 0.40 was used. 

^Totai-DB> A^si-DB> and Af c _ DB for several treatment tem- 
peratures are shown in Figure 3. Post-HPT, Si-QDSL de- 
fect density (1.1 x 10 19 cm" 3 ) clearly reduced compared 
with the defect density before HPT. The defect density 
for 200°C treatment is still high because hydrogen diffu- 
sion is insufficient. Hydrogen intrusion depth for 60-min 
HPT can be estimated to be below 100 nm, and a several 
dangling bonds remain in the deep area of the film. The 
defect density for 300°C treatment is lower than that at 
200°C. A large amount of hydrogen reaches the interface 
of the film and substrate during the 60-min HPT. The 
measured g value in this sample was 2.00275, which is 
quite similar to the g value of C-DB, meaning that 
Afei-DB is less than A/c-db- Based on Equation 5, A/si-DB 
is estimated to be 2.2 x 10 16 cm" 3 , indicating that 
Si-DBs can be efficiently passivated by the incorporated 



Table 1 Thicknesses estimated by fitting of the 
spectroscopic ellipsometry measurements of Si-QDSLs 



Parameters 


300°C 


400°C 


500°C 


600°C 


MSE 


11.56 


12.22 


13.37 


13.30 


7" s (nm) 


33.1 


11.5 


15.2 


6.5 


T (nm) 


167.7 


212.8 


224.7 


246.1 



hydrogen. For the 400° C treatment, defect density de- 
creases to 3.7 x 10 17 cm" 3 , which is comparable with the 
defect density of an a-SiC film. The g value for 400°C 
treatment was higher than that for 300°C treatment, 
indicating that C-DBs - which are dominant in the total- 
DBs - significantly decrease despite the increment in 
Si-DBs. For the 500°C treatment, defect density increases 
despite efficient hydrogen incorporation in the Si-QDSL, 
showing that the hydrogen atoms are thermally activated 
from the Si-H bond state to the interstitial state above 
300°C and from the C-H bond state to the interstitial 
state above 400° C. These temperatures mostly corres- 
pond to the temperatures of dehydrogenation from Si-H 
bonds and C-H bonds, which are approximately above 
300°C [26] and 450°C to 550°C [27], respectively. In the 
500°C treatment sample, a large amount of hydrogen 
atoms were in the interstitial sites; they did not contrib- 
ute to the passivation of the dangling bonds. 

Figure 4 shows the Raman spectra of the Si-QDSLs 
after 60-min HPT at different temperatures. The peak 
found between 2,000 and 2,100 cm" 1 corresponds to the 
Raman shift originating from the stretching mode of 
Si-H w bonds. The intensity of the peak from Si-H w bonds 
gradually weakens as the treatment temperature in- 
creases, indicating that the Si-H w bonds decomposed by 
the thermal activation of terminal hydrogen atoms. This 
trend agrees with the increment of Nsi-DB- 

During HPT, the surface is damaged by atomic hydro- 
gen. The thickness of the surface damaged layer is 
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Figure 6 RMS roughness measured by AFM and thicknesses of 
surface damaged layers estimated by spectroscopic ellipsometry. 
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Figure 7 Cross-sectional TEM images. At the near-surface of (a) 350°C treatment sample, (b) 600°C treatment sample, (c) magnified image of 
350°C treatment sample, and (d) magnified image of 600°C treatment sample. 



dependent on the treatment temperature. The thickness 
of the surface damaged layer was estimated by spectro- 
scopic ellipsometry. A schematic of the structure used 
for the analysis is shown in Figure 5. The TauoLorentz 
model was applied to the optical modeling of the Si- 
QDSL layer, and the surface damaged layer was assumed 
to be the effective medium approximation (EMA) layer 
in which 50% void exists. The estimated thicknesses of 
the Si-QDSL layers T, the thicknesses of the surface 
damaged layers T SJ and the mean square error (MSE) of 
each fitting are summarized in Table 1. T s of an as- 
annealed Si-QDSL was approximately 2 nm, while the 
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T s of the treated Si-QDSLs drastically increased, indicat- 
ing that the Si-QDSL structure in the surface region was 
broken by the atomic hydrogen. 

The thicknesses T and T s strongly depend on the 
treatment temperature. T decreases as the treatment 
temperature increases; this tendency is related to the 
hydrogen concentration at the near-surface for each 
treatment temperature. A large amount of hydrogen 
introduced into amorphous silicon contributes to the 
structural reconstruction by breaking the weak Si-Si 
bonds [28,29]. Further, surface morphologies were mea- 
sured by AFM. The root mean square (RMS) surface 
roughness of the samples is shown in Figure 6. RMS sur- 
face roughness is almost independent of the treatment 
temperature, whereas the damaged layer thickness 
measured by spectroscopic ellipsometry decreased with 
treatment temperature, indicating that HPT at low 
temperature introduces a damaged layer with lower 
refractive index than that of Si-QDSL. To investigate 
further, TEM observations of the Si-QDSLs were con- 
ducted. Figure 7a,b shows TEM images of the 350°C and 

Table 2 Thicknesses estimated by fitting of the 
spectroscopic ellipsometry measurements of 
surface-etched Si-QDSLs 



300 400 500 600 


Parameters 


300°C 


400°C 
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Temperature (°C) 
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14.72 


15.90 


Figure 8 RMS roughness measured by AFM and thicknesses of 
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the surface damaged layers of Si-QDSLs after RIE. 
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600°C treatment samples, and Figure 7c,d shows the 
magnified images of each sample. In the magnified 
images, existence of the Si-QDs is indicated using red 
circles. The irradiated electrons are transmitted through 
the sample without scattering in the white region, show- 
ing that the material density at the near surface is ex- 
tremely low in the white region. Detailed analysis of the 
TEM images revealed that the two periods of superlat- 
tice layers were completely removed by 350°C HPT. 
Two or three periods of superlattice layers were found 
to be damaged. On the other hand, for the 600°C treat- 
ment sample, no removal of the layers was observed 
during the HPT treatment; only the one-period superlat- 
tice layer was damaged. This result agrees with the 
thickness of the damaged layer estimated by the spectro- 
scopic ellipsometry. 

The damaged layer is defective and no longer acts as a 
Si-QDSL. Therefore, the existence of the damaged layer 
is a cause of the degradation of Si-QDSL solar cell per- 
formance. The removal of the damaged layer without 
additional damage is very important. Therefore, etching 
of the damaged layer was performed using RIE. RMS 
roughness measured by AFM and the damaged layer 
thicknesses estimated by spectroscopic ellipsometry of 
the Si-QDSLs after RIE are shown in Figure 8. The esti- 
mated thicknesses of the Si-QDSL layers T, the thick- 
nesses of the surface damaged layers T s , and the MSE of 
each fitting are summarized in Table 2. The observed 
RMS roughness was less than 3 nm, which was almost 
the same as that of the sample before RIE. The thick- 
nesses of the surface damaged layers estimated by spec- 
troscopic ellipsometry were almost the same as those of 
the RMS roughness. In general, surface roughness is also 
modeled using the EMA model for ellipsometry analysis; 
thus, the estimated T s reflects surface roughness, and no 
damaged layer exists on the surface. These results clearly 
indicate that RIE can remove the damaged layer without 
additional damage to the sample; RIE is therefore the 
key to improve the film quality of Si-QDSLs and the p/i 
interface in Si-QDSL solar cells. 

Conclusions 

Hydrogen plasma treatment temperature dependences 
of defect densities and hydrogen concentrations in 
Si-QDSLs as well as the surface morphologies of 
Si-QDSLs were investigated. Hydrogen could be quickly 
incorporated as the treatment temperature increases. 
On the other hand, dehydrogenation of hydrogen atoms 
terminating the dangling bonds is dominant during 
high-temperature treatments. The optimal treatment 
temperature was found to be approximately 400°C, and 
a defect density of 3.7 x 10 17 cm" 3 was achieved, which 
is comparable to the defect density of a typical a-SiC:H 
film. In addition, damaged layer was found to form on 



the surface by HPT; this damaged layer can be easily re- 
moved by RIE without additional damage to the sample. 
Thus, HPT and damaged layer removal process are very 
important for the fabrication of Si-QDSL solar cells. 
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